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ABSTRACT 



We show that the coronal heating and the acceleration of 
the fast solar wind in the coronal holes are natural conse- 
quence of the footpoint fluctuations of the magnetic fields 
at the photosphere by one-dimensional, time-dependent, 
and nonlinear magnetohydrodynamical simulation with 
radiative cooling and thermal conduction. We impose 
low-frequency (< 0.05Hz) transverse photospheric mo- 
tions, corresponding to the granulations, with velocity 
(dv±) = 0.7km/s. In spite of the attenuation in the 
chromosphere by the reflection, the sufficient energy of 
the generated outgoing Alfven waves transmit into the 
corona to heat and accelerate of the plasma by nonlin- 
ear dissipation. Our result clearly shows that the initial 
cool (10 4 K) and static atmosphere is naturally heated up 
to 10 6 K and accelerated to ~ 800km/s, and explain re- 
cent SoHO observations and Interplanetary Scintillation 
measurements. 

Key words: magnetic fields - plasma - magnetohydrody- 
namics - Sun : corona - solar wind - waves. 



1. INTRODUCTION 



It is still poorly understood how the coronal heating and 
the acceleration of the high-speed solar winds are accom- 
plished in the open coronal holes which mainly exist in 
the polar regions except at the solar maximum phase. The 
origin of the energy to heat and accelerate the plasma is 
believed to be in the surface convection. This energy is 
lifted up through the magnetic fields and its dissipation 
leads to the plasma heating. In general, the problem of 
the coronal heating and the solar wind acceleration is to 
solve how the solar atmosphere reacts to the footpoint 
motions of the magnetic fields. Then, an ideal way to 
answer the problem is to solve the transfers of mass, mo- 
mentum, and energy in a self-consistent manner from the 
photosphere to the interplanetary space, although such an 
attempt has not been successful yet. 
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In the coronal holes the Alfven wave is regarded to 
play an important role (e.g. Belcher 1971), since it 
can travel a long distance to contribute not only to the 
coronal heating but to the solar wind acceleration. The 
Alfven waves are excited by steady transverse motions 
of the field lines at the photosphere (e.g. Cranmer & 
van Ballegooijen 2005), while they can also be pro- 
duced b y contin ual reconnec tions above the photosphere 
(lAxford & Mckeinziel 119971) . The latter process might 
be responsible for generation of the high-frequency (up 
to 10 4 Hz) ioncycl otron wave highlighted for the heating 
of the heavy ions (iKohl et all fT998). However, it is diffi- 
cult to sufficiently heat the protons having the higher res- 
onant frequency by the high-frequency waves, because 
the wave energy is already absorbed by the heavy ions 
with hig her mass-t o-charge ratio (i.e. lower resonant fre- 
quency) (Cranmer, 2000). On the other hand, the low- 
frequency ( ~ 0.1Hz) Alfven wave does not have such 
a drawback; it can propagate a long distance and heat the 
bulk of the plasma. Since we focus on the heating of 
the major part of the plasma, we study the low-frequency 
Alfven waves by the steady footpoint fluctuations. 

The low-frequency Alfven waves have been intensively 
studied in the context of the heating a nd acceleration 
of the sola r wind plasma . For example, lOughton et alJ 
(1200 ll) and lDmitruk et alJ (120021) examined turbulent cas- 
cade o f the Alfven wa ves in tangential directions. Re- 
cently, lOfmanl (120041) performed three-fluid and two- 
dimensional simulations for the wave-driven fast solar 
wind. These works are quite important in understand- 
ing the physical mechanism of the solar wind accelera- 
tion. One of the shortcomings of the previous research is 
that the most calculations for the wave-driven solar wind 
adopt the fixed 'coronal base', instead of the photosphere, 
as the inner boundary. In reality, however, the corona dy- 
namically interacts with the chromosphere located below 
by chromosp heric evaporation ( Hammer, 1982) and wave 
transmission ( Kudoh & Shibata, 1999) so that the coronal 
base varies time -depende ntly. Alth ough some calcula- 
tions tLie-Svendsen et all 1200 It 120 02) include the chro- 
mosphere and the transition region to properly take into 
account these effects, they assume unspecified mechan- 
ical energy deposition which requires an ad hoc heating 
function with the phenomenological dissipation length. 
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Figure 1. Schematic picture of our simulation (especially 
for NON-solar physicists). We consider the waves in the 
open flux tube rooted at the photosphere in the coronal 
hole to 0.3 AU. Note that a coronal hole (dark region) is 
clearly seen in the north polar region in Yohkoh/SXT im- 
age. 



In contrast, we self-consistently treat the transfer of 
the mass/momentum/energy by dynamically solving the 
wave propagation from the photosphere to the interplan- 
etary region. We handle the heating and acceleration by 
the waves in an automatic way without the phenomeno- 
logical heating function. Our aim is to answer the prob- 
lem of the heating and ac celeration in the coronal holes 
by the forward approach (Gudiksen & Nordlund, 2005) 
with the minimal unknown parameters. 



2. SIMULATION METHOD 



We consider one-dimensional open flux tube which is 
super-radially open (fig[D, measured by heliocentric dis- 
tance, r. The simulation region is from the photosphere 
(r = li?s) with density, p = 10 _7 g cm -3 , to 65i?s 
(0.3AU), where i?s is solar radius. Radial field strength, 
B r , is given by conservation of magnetic flux as 



B r r f — const. 



(1) 



where / is a super-radial expansion factor 
( Kopd & O rraljl. Il97fil) . In this paper B r is set to 
be 161G at the photosphere, 5G at low coronal height, 
r = 1.02i? s , and B r = 2.lG(R s /r) 2 in r > 1.5R S 
by controlling /. We give transverse fluctuations of the 
field line by the granulations at the photosphere, which 
excite Alfven waves. We consider the fluctuations with 



6 x 10~ 4 < v < 0.05Hz (period of 20seconds — 30min- 
utes), and root mean squared (rms) average amplitude 
(dv±) ~ 0.7km/s corre sponding to observed velocity 
amplitude ~ lkm/s (Holwege r. Gehlsen. & Rulandl 
1978). At the outer boundary, non- reflecting co ndition 
is imposed for all the MHD waves (Suzuki & Inutsuka, 
2005a b), which enables us to carry out simulations for a 
long time until quasi-steady state behaviors are achieved 
without unphysical wave reflection. 

We dynamically treat propagation and dissipation of the 
waves and heating and acceleration of the plasma by solv- 
ing ideal MHD equations with the relevant physical pro- 
cesses : 
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power spectrum, P(v) 



in frequency between 



rV/ dr 

where p, v, p, e, B are density, velocity, pressure, 
specific energy, and magnetic field strength, respec- 
tively, and subscript r and _L denote radial and tangen- 
tial components. G and Ms are the gravitational con- 
stant and the solar mass. F r is thermal conductive flux 
and or is radiative cooling jlLandini & Monsignori-Fos'sil 
1990; lAnderson & Athavl Il989t iMorivasu et all 12004^ 
We adopt 2nd-order MHD-Godunov-MOCCT scheme 
for updating the physical quantities (Sano & Inutsuka 
2005), of which an advantage is that no artificial vis- 
cosity is required even for strong MHD shocks. We fix 
14000 grid points with variable sizes in a way to resolve 
Alfven waves with v = 0.05Hz by at least 10 grids per 
wavelength to reduce numerical damping. 

The advantages of our simulation are (i) automatically 
treating wave propagation/dissipation and plasma heat- 
ing/acceleration without an ad hoc heating function (ii) 
with minimal assumption, i.e., the transverse footpoint 
motions in the given flux tube (iii) in the broadest re- 
gions with respect to the density contrast amounting to 
15 orders of magnitude from the photosphere to 0.3AU. 
On the other hand, the disadvantages are due to the 1- 
dimensional MHD approximation; we assume that the 
wave propagation is restricted to the direction parallel to 
the field line and the plasma behaves as the one-fluid with 
the Boltzmann particle distribution. We will discuss this 
issues later. We would like to emphasize however, that 
this is the most self-consistent simulation for the solar 
wind acceleration, in spite of these shortcomings. 
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Figure 2. Time evolution of the atmosphere. Outflow 
speed, v r (km/s) (upper-left), temperature, T(K) (upper- 
right), density, rms amplitude of transverse velocity, 
(dv±) (km/s) (lower-left), and p(g/cm 3 ) (lower-right) are 
plotted. Black dashed, green, blue, and red solid lines 
are results at t = 0, 20, 340, & 2573 mins., respectively, 
whereas results at t — do not appear in (dv±) because 
it equals to 0. The results are averaged by integrating for 
3 minutes to take into account effects of the exposure time 
for comparison with observations. 



3. RESULTS 



We initially set static and cool atmosphere with temper- 
ature, T = 10 4 K; we do not impose the corona and the 
solar wind at the beginning. Figure|2j shows how the 
coronal heating and the solar wind acceleration are real- 
ized by the generated low-frequency Alfven waves. We 
plot v r (km/s), T(K), p(g/cm 3 ), and (dv±)(km/s) aver- 
aged from v± as a function of (r — Rs)/Rs at differ- 
ent time, t = 0, 20, 340 & 2573 minutes. As time goes 
on, the atmosphere is heated and accelerated effectively 
by dissipation of the Alfven waves. Temperature rises 
rapidly in the inner region even at t — 20minutes, and the 
outer region is eventually heated up by both outward ther- 
mal conduction and wave dissipation. Once the plasma is 
heated up to the coronal temperature, mass is supplied to 
the corona mainly by chromospheric evaporation due to 
the downward thermal conduction. This is seen in tem- 
perature structure as an inward shift of the transition re- 
gion, which is finally located around r = 6 x 10 _3 i?s 
(~4000km), whereas it moves up and down by the time- 
dependent chromospheric evaporation (heating) and the 
wave transmission. As a result, the coronal density in- 
creases by two orders of magnitude. While the wind ve- 
locity exceeds lOOOkm/s at t = 340minutes on account 
of the initial low density, it gradually settles down to 



< lOOOkm/s as the density increases. (dv±) also settles 
down to the reasonable value at the final stage. Temper- 
ature structure is smoother owing to the thermal conduc- 
tion than the other quantities showing fluctuated behav- 
iors due to the waves. 

We have found that the plasma is steadily heated up 
to 10 6 K in the corona and flows out as transonic wind 
with v r ~ 800km/s at the outer boundary (=0.3AU) 
when the quasi steady-state behaviors are achieved after 

t ~ 1800minutes. This is the first numerical simu- 
lation which directly shows that the heated plasma actu- 
ally flows out as the transonic wind, initiated from the 
static and cool atmosphere, by the Alfven waves. The 
outflow speed becomes up to lOkm/s even below 10 4 km 
above the photosphere (c.f. Tu et al.2005), although it 
is difficult to distinguish from the longitudinal wave mo- 
tions. The sonic point where v r exceeds the local sound 
speed is located at r ~ 2.5i?s and the Alfven point 
is at r ~ 24i?s. Obtained proton flux at 0.3AU is 



(2 ± 0.5) x 10 (cm s ), corresponding to 



N p v ~ (1.8 ± 0.5) x lO^cm-V 1 ) at 1AU for N p v oc 
r~ 2 , which is consistent with the observed high-speed 
stream around the earth (Asc hwanden. Poland. & Robin, 
2001). 

In fig|3]we compare the result at t = 2573minutes with 
recent observation in the high-speed solar winds from 
the polar regio ns by Sola r & He li ospheric Obs er vatory 
fSoHO) (IZangrilli et all l2002t [Teriaca et all l2003| 
Fludra, Del Zanna, & Bromage, 1999; Lamv et al„ 
19971 IWilhelm et all 1 19981 iBaneriee et all I199E 



'Movie file is available : 
http://www-tap.scphys.kyoto-u.ac.jprstakeru/research/suzuki_2005()6.mpg 



Esser et all 1 19991) and Interplanet a ry Scintil lation (IPS) 
measuremen t s (iGrall et all [1996; Koiima et all 120041 
ICanals et all 120021) . The figure shows that our forward 
simulation naturally form the corona and the high-speed 
solar wind which are observed except small differences in 
detailed structures. Observed outflow speed in the inner 

corona(< 3i?s) and outer region( ~ 20i?s) are mostly 
explained by our simulation within the obse rved errors . 
Some of the observed data around r ~ lOfk JGrall et all 
1996) exceed our result, whereas it is reported that 
these data migh t reflect wave phenomen a rather than the 
outflow speed ( lHarmon & Coles 1 120051) . The simulated 
temperature shows a decent agreement with the elec- 
tron temperature in the inner region by SoHO (Fludra 
et al.1997). Density and transverse amplitude show 
reasonable agreements with the observations. 

Our result manifestly shows that the heating and acceler- 
ation of the high-speed winds from the open field regions 
can be almost completely explained by the MHD dissi- 
pation mechanisms of the low-frequency Alfven waves. 
The result is quite convincing since we automatically 
solve the transfers of mass/momentum/energy and the 
propagation of the Alfven waves excited at the photo- 
sphere without any ad hoc assumptions. 

Figure^ presents energy flux of outgoing Alfven wave, 
incoming Alfven wave, and outgoing MHD slow (sound) 
wave at t = 2573minutes. To focus on the amount of dis- 
sipation the energy flux is normalized by the cross section 
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Figure 3. Comparison of the result at t — 2573mins (red solid lines) with observations summarized below. Quantities in 
four panels are the same as in fig^\except lower-right panel showing electron density, N e (cm~ 3 ), instead of p(g cm~ 3 ). 
Scale s in both horizontal and vertical axises are changed from fig]2\ a: Vertical error bars with green x's \Teriaca et all 
2003) and blue triangles iZans rllli et all \2002\) are proton outflow speed in polar regio ns by SoHO. Blue sq uare with 
errors is velocity by IPS measurements averaged in 0.13 - 0.3AU of high-la titude regions\Konma et~al\ ^2004). Crossed 
bars with and without circles are measurements by VLBI and IPS( EISCAT) ( Grail et al., 1 996 ). Vertical error bars with 
circles are data based on observation by SPA RTAN 201-01 (Habbal et all \l 99?\) . b: Pink circles are electron temperature 
by CDS/SoHO Wludra. Del Zanna. & Brom aee\\19 9~<kl. c: Pink circles and blue stars are obser vations by SUM ER/SoHO 
\Wilhelm e^d^ and b y CDS/SoHO \Teriaca et all 12003). respectively. Green triangles \Teriaca et all \2003) and 

squares $Lamvetall\199'A) are obs ervations by LASCO/S oHO. d: Blue circles are non-thermal broadening inferred 
from SUMER/SoHO measurements \Baneriee et all \199li). Cros s hatched region is empirically derived non-thermal 
broadening based on UVCS/SoHO observation iEsser et qf\\199 9 \ ). Green er ror bars with circles are transverse velocity 
fluctuations derived from IPS measurements by EISCAT ( Cana lset all \2Q02). 
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Figure 4. Energy flux of outgoing Alfven (red solid), in- 
coming Alfven (blue dotted), and outgoing MHD slow 
(green dashed) waves at t — 2573mins. Hatched region 
indicates the chromosphere. The energy flux is normal- 
ized by the cross section (A = r 2 f) of the flux tube at 
r = 1.02i?g and A c denotes A at that location. 

of the flux tube at r = 1.02i?g where it expands super- 
radially by 30 times from the photosphere. Note that 
the real energy flux in the lower (upper) regions is larger 
(smaller). The outgoing and incoming Alfven waves are 
decomposed by correlation between v± and B±_ (Elsasser 
variables). Extraction of the slow wave is also from fluc- 
tuating components of v r and p. The derived energy flux 
is averaged by integrating for 30 minutes to smooth out 
variation due to phase. 

The figure shows that the outgoing Alfven wave dissi- 
pate quite effectively; less than 10 -3 of the initial en- 
ergy is remained as that associated with the waves at the 
outer boundary. First, a sizable amount is reflected back 
downward below the coronal base (r — i?g < 0.01-Rg), 
which is clearly illustrated in the energy flux of the in- 
coming Alfven wave following that of the outgoing com- 
ponent with slightly smaller level. This is because the 
wave shape is considerably deformed owing to the steep 
density gradient; a typical variation scale (< 10 5 km) of 
the Alfven speed becomes comparable or even shorter 
than the wavelength (= 10 4 - 10 6 km). Although the 
energy flux, ~ 5 x 10 5 erg cm _2 s -1 , of the outgoing 
Alfven waves which penetrates into the corona is only 
~ 15% of the input value, it meets the requirement for the 
energ y budget in the coronal holes ( Withbro^^^oyei, 

Second, MHD slow (sound) waves (Sak uraietall f2002) 
are generated in the corona as shown in the fig- 
ure. The amplitude of the Alfven waves is amplified 
through the upward propagation, as a result, a non- 
linear term due to variation of magnetic pressure ac- 
companying wi th the Alfven waves ex cites longitudi- 
nal slow waves(Kudoh & Shibata, 199S>b- They eventu- 
ally steepen to form MHD slow shocks llSuzukil 120021) 



which efficiently convert kinetic and magnetic energy 
to heat. The coronal heating and wind acceleration are 
thus far achieved by transfer from the energy and mo- 
mentum of the outgoing Alfven waves. Linearly polar- 
ized Alfven waves directly steepen in themselves to form 
MHD fast shocks faollweel fl982l ISuzukil l200l . which 
also contributes to the heating, though it is less domi- 
nant. The incoming Alfven waves are generated in the 
corona by the reflection of the outgoing ones by the den- 
sity fluctuations due to the slow waves. The reflected 
waves further play a role in the dissipation of the out- 
going Alfven waves by nonlinear wave-wave interaction. 



4. SUMMARY AND DISCUSSIONS 



We have performed one-dimensional MHD simulation 
for the low-frequency Alfven waves excited at the pho- 
tosphere. We incorporate radiative cooling and thermal 
conduction, and self-consistently treat the mass, momen- 
tum, and energy transfer by nonlinearly and dynamically 
solving the wave propagation and dissipation. The ad- 
vantage of our simulation is the forward approach with 
the minimal assumptions. 

Our result has manifestly shown that the formation of the 
corona and the acceleration of the fast solar wind in the 
coronal holes are the natural outcome of the footpoint 
motions of the field line. The photospheric fluctuations 
generate the low-frequency Alfven waves which propa- 
gate upwardly. About 15% of the initial Alfven waves 
in energy flux can transmit into the corona after surviv- 
ing the non-WKB reflection in the chromosphere and 
the transition region. The amplitudes of these out-going 
waves are amplified due to the density stratification so 
that they effectively dissipate by the nonlinear effects and 
heat and accelerate the coronal plasma. 

A key mechanism in the coronal heating and the solar 
wind acceleration is the generation of the MHD slow 
wave. Thus, one of the predictions from our simulation is 
the existence of the longitudinal fluctuations in the solar 
wind plasma. This is directly testable by future missions, 
Solar Orbiter and Solar Probe, which will approach to 
~45 and 4 i?g, respectively, being in our simulation re- 
gion. They can determine how much fractions of the fluc- 
tuations are in the transverse and longitudinal modes by 
in situ measurements which are directly compared with 
our result. 

In this paper we have considered the wave prop- 
agation and dissipation in one-dimensional MHD 
approximation. However, its validity needs to be 
examined by future studies. It is important to in- 
vestigate the m ultidimensional effects for the waves 
(OfmaiJ l200l . such as refraction (Boadan et al., 
2003), phase mixing! Hevvaerts & Priest, 1983; 
Nakariakov, Roberts, & Mura wskit 1 19981). turbulent 
casca de into the transvers e direction Joughton et all 
1200 It iDmitruk et all [2002). The waves could suffer 
collisionless damping (e.g. Suzuki et al.2005), and 
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this mechanism might also modify the dissipation rate. 
Kinetic effects of the multicomponent plasma might be 
also important especially in the outer region with low 
density. 

We need to check how much the wave dissipation and 
the consequent plasma heating are modified by these pro- 
cesses. However, we cannot study them in the global sim- 
ulations as performed here because of the huge density 
difference. Hence, we should firstly investigate the de- 
tailed processes in the local simulations, and then, com- 
pare with global results; the global simulation in the one- 
dimension and the local simulations in the two- or three- 
dimensions should play complimentary roles. 
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